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Abstract 

a^"■ 
a^ : 

' We perform a first principles computation of the spectrum of gravity waves 

. produced in open inflationary universes. The background spacetime is taken 

! to be the continuation of an instanton saddle point of the Euclidean no bound- 

^ I ary path integral. The two-point tensor correlator is computed directly from 

^f) ' the path integral and is shown to be unique and well behaved in the in- 

frared. We discuss the tensor contribution to the cosmic microwave back- 
^ . ground anisotropy and show how it may provide an observational discriminant 

, between different types of primordial instantons. 

o ■ 
cn ■ 

O ■ I- INTRODUCTION 

o\ ' 

. The inflationary universe scenario provides an appealing explanation for the smoothness 

p ^! and flatness of the present universe, as well as a mechanism for the origin of density fluc- 
^ I tuations. Until recently it was believed that inflation inevitably predicted a flat Qq = 1 
i-^ I universe. However in ]I| it was shown that with mild fine tuning an open universe is also 
^ ■ possible. The potential must have a sharp false vacuum in which the field is assumed to 
^ . have become trapped. The field is then assumed to tunnel out via an instanton known as 
• 1— I ! the Coleman-De Luccia instanton 0, producing a bubble within which slow roll inflation 
^ I occurs. The interior of the bubble produced via the Coleman-De Luccia instanton is an 
I infinite inflating open universe. 

Such models provide important counter-examples to the standard folklore but require 
quite contrived scalar field potentials. Recently, however. Hawking and one of us showed that 
open inflation can occur much more generally. We found a new class of instantons that 
exist for essentially any inflationary potential, and provide saddle points of the Euclidean 
path integral. The continuation of these instantons is similar to that of the Coleman- 
De Luccia instantons, and they define initial conditions for open inflationary universes. 
Although the Hawking- Turok instantons are singular, the singularity is mild enough for the 
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quantization of perturbations to be well posed 0, fH, 0- In this paper we compute the 
spectrum of gravity waves for both Coleman-De Luccia and Hawking- Turok instantons. 

This paper is a companion to ref. 0, where scalar fluctuations about open inflationary 
instantons were calculated. Here we perform the analogous calculation for the tensor fluctu- 
ations and discuss possible observable signals in the CMB anisotropy power spectrum. The 
calculation is performed in the framework of the Euclidean no boundary proposal due to 
Hartle and Hawking as discussed in [^]. The correlator is computed in the Euclidean re- 
gion where it is uniquely determined by a Gaussian integral, and then analytically continued 
in the coordinates of the classical background solution into the Lorentzian region of interest. 
Our main result for the tensor correlator (^) is given in a form which is straightforward 
to compute numerically. We defer detailed numerical calculations of the CMB anisotropies 
to a future paper ^ in which both scalar and tensor contributions for a variety of scalar 
potentials will be discussed. 

There have over the last few years been many papers exploring similar calculations, 
mostly making one approximation or another |10|, [Tl|. Very recently, Garriga et. al. 



have independently obtained formulae for the scalar and tensor correlators similar to ours 
1^ . These formulae have been numerically implemented in ref. which gives results for 



some examples of Coleman-De Luccia instantons calculated without approximation. 

We feel that the derivation given here is significantly clearer than in these papers, and 
that our method has several important conceptual advantages. First, all earlier authors 
performed a mode by mode analysis. In this framework, one requires a prescription for the 
vacuum state for each perturbation mode and this is taken to be the state in which the 
positive frequency modes are regular on the lower half of the instanton. This prescription 
is rather ad hoc. In contrast, our method is to simply perform the Euclidean no boundary 
path integral. This automatically gives a unique Greens function. There is no need for an 
additional prescription, indeed imposing one is contrary to the spirit of the no boundary 
proposal (see the discussion in |^). The whole idea of the Euclidean no boundary proposal 
is that an essentially topological prescription should define the initial state of the universe. 
Analyticity arises because the background solution is a solution to a differential equation. 
Divergent fluctuation modes have infinite Euclidean action and are therefore suppressed 
in the path integral. Second, in the matching method of Garriga et. al., they devote a 
great deal of effort to determining the action for perturbations in region II, the part of the 
Lorentzian spacetime exterior to the open universe region. This introduces considerable 
technical complexity since the spatial hypersurfaces used in their canonical quantisation 
approach are inhomogeneous in region II. Our approach is to analytically continue directly 
from the Euclidean region into the open universe. Region II is just a part of the continuation 
route with no special significance. Third, as emphasised in 0, we deal throughout directly 
with the real space correlator. In this approach 'super-curvature' modes are automatically 
included and their relation to the 'sub-curvature' modes is thereby made clear. A related 
fact is that we find that the real space correlator to be infrared finite even in perfect de 
Sitter space, as mentioned below. Finally, Garriga et. al. only give formulae for equal-time 
correlators. To compute the microwave anisotropies one requires the unequal-time correlator, 
which we give here. We are also careful to define the continuation of the conformal time 



coordinate into the Euclidean region, which is not explained in ref. [14 



The paper is organised as follows. In section 2 we describe the relevant path integral 
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and the model-dependent Schrodinger operator which occurs in the Euchdean action. We 
show that for singular instantons the singularity acts as a reflecting boundary, flxing Dirichlet 
boundary conditions for the perturbation modes Q], The Euclidean tensor correlator 
is computed from the path integral in section 3. In this calculation we need several properties 
of maximally symmetric bitensors on 5*^, which are described in Appendix A. Section 4 
describes the analytic continuation to the open universe. Finally, section 5 is devoted to 
the Sachs- Wolfe integral to determine the contribution of gravitational waves to the CMB 
anisotropy. Here we comment on possible observational distinctions between Coleman-De 
Luccia and Hawking- Turok instantons. 

We conclude this introduction with two technical remarks. First, the question of discrete 
'supercurvature' modes arises in the tensor calculation just as in the scalar case p. Here, 
however, we find that although the relevant Schrodinger operator possesses a bound state 
just as in the case of scalar perturbations here it does not generate a 'super-curvature 
mode'. Instead the relevant mode is a time-independent shift in the metric perturbation 
which may be gauged away. This is in agreement with refs. |T^, |12|. Second, it has been 



claimed in the much of the previous literature that the spectrum of gravity waves in pure 
de Sitter space is infrared divergent p but that the divergence disappears once the 



existence of the bubble wall is taken in account [|TU[. In our approach we find a different 
result. Neglecting the gauge mode previously mentioned, the two point correlator has a well 
defined long-wavelength hmit even in perfect de Sitter space. We shall investigate this issue 
further in future work. 



II. THE PATH INTEGRAL FOR TENSOR FLUCTUATIONS 

In quantum cosmology the basic object is the wavefunctional \E' [hij,(f)], the amplitude 
for a three-geometry with metric hij and field configuration 0. It is formally given by a path 
integral 

^[h„<P]^ f'"''[Dg]me''^^''^\ (1) 

Following Hartle and Hawking 0] the lower limit of the path integral is defined by continuing 
to Euclidean time and integrating over all compact Riemannian metrics g and field config- 
urations (p. If one can find a saddle point of (|I]), namely a classical solution satisfying the 
Euclidean no boundary condition, one can in principle at least compute the path integral as 
a perturbative expansion to any desired power in h. 

In this paper, we shall compute the two-point tensor fiuctuation correlator about classical 
solutions describing the beginning of open infiationary universes, to first order in h. The 
principles are described in namely that we compute the correlator in the Euclidean 
region where the exponent iS in the path integral becomes —Se = —{So + S2), where Se is 
the Euclidean action. So is the instanton action and S2 the action for fiuctuations. We shall 
keep the latter only to second order, this being all that is needed to compute the quantum 
fiuctuations to leading order in h. The correlator is then given by a Gaussian path integral 

{U,[x)U,j,[x ))- ^ ^^^^^ ^^^^^ . [Z) 
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The Lorentzian correlator is then obtained by analytically continuing in the coordinates of 
the background classical solution, into the open inflating region. 

The 0(4) symmmetric instantons of interest possess a line element of the form 
da^ + h'^{a)dVL\ where dVt\ is the line element on S^. Both Hawking- Turok and Coleman-De 
Luccia instantons possess a regular pole which we take to be at a = 0. As a approaches zero, 
we have h{a) a. The Coleman-De Luccia instantons have a second regular pole where 
h (Jm — o" where is the maximum value of a. In contrast Hawking- Turok instantons 
have h —>■ {a^n — cr)3 as ex — *• am- It is useful in both cases to introduce a conformal spatial 
coordinate satisfying dX = da/b{a), so that the line element takes the form 

ds^ = b\X) (dX^ + dnf) (3) 

For Hawking- Turok instantons we define 

X= / 7T^. (4) 

SO X = corresponds to the singular pole and X oo to the regular pole. For Coleman-De 
Luccia instantons X may be conveniently defined by ^^^^^ value of sigma 

for which is a maximum, and then X ranges from —oo to -|-cxo. We write the perturbed 
line element and the scalar field as 



ds'^ = b\X) ((1 + 2A)dX^ + Sidx'dX + {jij + hij)dx'dx^] 
= 0o(X) + 50. (5) 



and decompose Si and hij as follows WE 



Si = B\, + Vi. (6) 

Here A3 is the Laplacian and \ j the covariant derivative on the three-sphere. With respect to 
reparametrisations of the three-sphere, h, B and E are scalars, Vi and Fi are divergenceless 
vectors and tij is a transverse traceless symmetric tensor. 

One may expand the spatial part of each of these spin-r fields in terms of a complete set 
of harmonics, labelled by the eigenvalues of the Laplacian on 5^, Ap = + (r + 1) where 
p = in and n is an integer. In general the decomposition of a metric perturbation hij into 
a scalar, vector and tensor part is unique. Hence one can write E, Fi and tij back in terms 
of hij |T6[. For scalar = —4 and vector = —4 harmonics however, the decomposition 
is not unique and there appears a degeneracy between scalar- or vector-type perturbations 
and Pt=0 and pf = —1 tensor modes respectively. Treatment of the former is complicated 
by the involvement of the scalar field ^ , but the latter mode is unambiguously pure gauge. 
We will return to this point in section 4. 

The Euclidean action is 

S=^J d'^Va [-R + ^ V^0V> + r (0)) - i 1 d'x,/^K, (7) 
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where the surface term is needed to remove second derivatives. Substituting the decomposi- 
tion (P) into the action (|^), we keep all terms to second order. The scalar, vector and tensor 
quantities decouple. The scalar perturbations are studied in ref. 0. The vector perturba- 
tions are uninteresting to first order in % since they are forced to be zero by the Einstein 
constraints. The tensor perturbations give the following second order positive Euclidean 
action 



(8) 



where prime denotes differentiation with respect to the conformal coordinate X. If one 
performs the rescaling tij = b{X)tij, and integrates by parts one obtains 



S2 = ^l d^x^U, {k + 3- A3) - ^ / d'x^U,f^^-{X) 



(9) 



where the Schrodinger operator 



K 



b^ 

"dX^ ^ J 



1 = 



dX^ 



+ U{X). 



(10) 



The form of the potential U{X) is shown in Figure 1 for de Sitter space, as well as 
examples of a Coleman-De Luccia instanton and a Hawking- Turok instanton. The operator 
K has in all three cases a positive continuum starting at eigenvalue — 0, as well as a 
single bound state tij = b{X)qij{Q) at p = i. 

For singular instantons the surface terms in (|^) play a crucial role. The potential U {X) — >■ 
~IW as X — s> 0. The eigenmodes of K behave as X^/^ or X^/^lnX near the singularity. 
The latter modes contribute positive infinity to the surface terms in (j^), and are therefore 
suppressed in the path integral. Hence we see that as in the scalar case 0, the path 
integral unambiguously specifies the allowed fluctuation modes as those which vanish at the 
singularity. 



III. THE EUCLIDEAN GREEN FUNCTION 

To evaluate the path integral (0), we first look for the Green function G^"'*-' of 
the operator in (P). The Euclidean fluctuation correlator (|^) will then be given by 
b^^{X)b~^{X')G^-'^ . The Euclidean Green function satisfies 

^ (i: + 3 - A3) G^^,,,(x, X', n') = 6{x - x')^~U'\,^,{n - n'). (ii) 

If we think of the scalar product as defined by integration over 5*^ and summation over 
tensor indices, then the right hand side is the normalised projection operator onto transverse 
traceless tensors on S^. Since the eigenmodes of the Laplacian form a complete basis, we 
can write the last term as 
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FIG. 1. Potential U{X) occurring in the Schrodinger operator governing tensor perturbations 
about the various instanton solutions discussed in the text. The dashed line shows the potential 
for an S"^ instanton corresponding to perfect de Sitter space, where U{X) = -2/cosh2(X). The 
upper solid line shows the potential for a Coleman-De Luccia instanton, where — oo < X < oo, and 
the lower solid line that for a Hawking- Turok instanton, with singularity indicated by the vertical 
dotted line. The potentials have been shifted in X so their minima coincide. All three are very 
similar to the right of the minimum. To the left, the Hawking- Turok potential diverges as one 
approaches the singularity. The potential is reflectionless in the S'^ case, weakly reflecting in the 
Coleman-De Luccia case and totally reflecting in the Hawking- Turok case. 



s\,A^ - fi') = E E E E (12) 

k 'P=e,o 1=2 m=—l 

where on we have 

^SQvim = ^kqvim (13) 

with Afc = —k{k + 2) + 2. Here V = {e, o} labels the parity, the angular momentum on 
takes the values k = 2,3,4, .. and 2 < I < k and —l<m<l are the usual quantum 
numbers on the two-sphere. Note that (/ > 2) because a spin-two field has no monopole or 
dipole components. The eigenmodes are normalized by the condition 

/ Vld'xqf,'Jlq^:}U^ = 6'^' 6rv'6iv6^m' (14) 

The set of eigenmodes form a representation of the symmetry group 5*0(4) of the man- 
ifold. It follows in particular that the sum over P, I and m defines a maximally symmetric 
bitensor IITB 



,M = E <^i^^''^m (15) 



vim 



that depends only on the geodesic distance 0!) between the two spacetime points. 
The Green function G^\iji can only be a function of fi') if it is to be invariant under 
isometrics of the three-sphere. Note that the indices i,j lie in the tangent space over the 
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point Q while the indices lie in the tangent space over the point Q'. The general form 
of the bitensor ly^^^^ -ij, appearing in tensor fluctuation correlators has been obtained by 
Allen |]I9[ and is given in Appendix A below. Here we note already that in terms of the 
label p = i{k + 1) = in, the bitensor on has precisely the same formal expression as the 
corresponding object on H"^. Since we would like to analytically continue our result for the 
Euclidean two-point correlator into the open universe, we will use the label p = in from now 
on. We now return to equation (|TTp for the the Euclidean Green function. 
By substituting the following ansatz for the Green function, 

X, X') = AkY: G,{X, X')W^;^ (16) 

p=in 

into (0) and noting that in terms of p = in, we have Xk = p"^ + 3, we obtain an equation 
for the model-dependent part of the Green function, 

(k - p^) G,{X, X') = 5{X - X') (17) 
Let us first discuss the case of singular instantons. The solution to equation ( [T7| ) is 

G,{X,X') = [^;{X)%{X')Q{X - X') + ^>-^{X)^;{X')Q{X' - X)] , (18) 

where is the solution to the Schrodinger equation that goes as X^/"^ as X — and 

(X) is the solution going as e*^"^ = e~^-^ as X tends to infinity. The factor Ap is the 
Wronskian — of the two solutions. 

We shall ultimately be interested in re-expressing this solution as an integral over real 
values of p in order to continue it to the open universe. To do so we must extend the 
solutions defined above at p = in into the complex p-plane. (X) becomes \l/p(X), 
defined for all complex p to be the solution which tends to X^ as X ^ 0. Being a solution 
of a regular differential equation this is analytic for finite p in the complex p-plane. On the 
other hand, (X) is the analytic continuation of gp{X), defined on the real p axis to be 
the solution tending to e^^^ as X — >■ oo. This is the Jost function, and is analytic in the 
upper half p-plane [^. The two solutions may be expressed in terms of each other as 

(X), (19) 

and their Wronskian Ap = "^'pQp — g'p^p = —2ipa_p, independent of X. This too is analytic 
in the upper half p-plane. Zeros of a_p in the upper half p-plane correspond to normalisable 
bound states. They can only occur on the imaginary p-axis, and in the case of interest here 
the only zero in the upper half p-plane is at p = i. This zero corresponds to the bound state 
mentioned above. For X > X' we have the Green function as a discrete sum 



+ioo 

Gr''if^,X,X') = 4/. E 7^^t(X)%(X')Wir'(^') (20) 

p=3i ^P"- 



For regular Coleman-De Luccia instantons a similar procedure may be followed. Here X 
ranges from — oo to +oo and we define the two linearly independent mode functions g^p^^ (X), 
which tends to e~*^"^ as X —oo, and (7"^^* (X), which tends to e*^"^ as X — > 00. These can 
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be shown to be orthogonal and analytic in the upper half p-plane. As X — +00, we have 
gf^ (X) ^ CpC'P^ + dpC-^P^. Hence, the Wronskian \ = gf''g'j^^' - g'j^^^'gp^' = -2ipdp 

and the Green function G^* (/i, X, X') may be expressed in a form analogous to that for 
singular instantons. 

Before proceeding to the analytic continuation, let us demonstrate that our Euclidean 
Green functions are regular at the regular pole. This is a nontrivial check because the 
coordinates a and X are singular there, and the rescaling becomes divergent too, b{X) ~ 
cr~^ ~ e~^^ . In the large X, X' limit, (pO|) becomes 

Gf'^'ift, X, X') = 2«: £ i fe-"(^-^') + ^e-"(^+^')) (/i) (21) 



n=3 



For n > 3 the Gaussian hypergeometric functions F(3 + n, 3 — n, 7/2,2) that constitute 
the bitensor have a series expansion that terminates, and they essentially reduce to 

Gegenbauer's polynomials 0^1^(1 — 2z). Using then the identity 



5:cr(x)g' = (l-2xg + g2)"' (22) 



1=0 



with q = e i^^^')^ one easily sees that the sum (21) indeed converges 



We have the Euclidean Green function defined as an infinite sum ^0|). We wish to 
represent it as an integral over p. To do so we must extend the summand into the upper half 
p-plane. We have already defined the wavefunctions for all complex p but we need to extend 
the bitensor as well. When the Green function is expressed as a discrete sum, it involves the 
bitensor (yu) evaluated at p = in with n integral. At these values of p, the bitensor is 

regular at both coincident and opposite points on S^, that is at yU = and /j, = tt. However, 
if we extend p into the complex plane we lose regularity at fi = 0. This is clear from (^Tj). 
For if we distort the p integral to run along the real axis, and use the completeness relation 
for the eigenfunctions ippi^X), it follows that W^(p)'*"' (/u) obeys a differential equation with a 
delta function source at = (see the discussion of the scalar case in p). Similarly, when 
we extend (//) into the complex p-plane, we must maintain regularity at /i = vr, since 

there is no delta function source there. 

The condition of regularity at vr imposed by the differential equation for the Green func- 
tion is sufficient to uniquely specify the analytic continuation of Wf^^^^ ■' (/i) into the complex 
p-plane. To see this, we note from Appendix A that The bitensor involves coefficent functions 
a and f3 which are hypergeometric functions of the variable z = cos^(/i/2). For coincident 
points, z = 1 but for antipodal points z = 0. There are two independent solutions of the 
hypergeometric equation, namely a{z) and a{l — z). They are related by the transformation 
formula (eq. [15.3.6] in 



7 7 
2^1(3 + 3 - ip, -, z) = (- coshj97r)2Fi(3 + ip,3- ip, z) 

r(i)r(f) , 5 ^,1 . 1 . 3 ^ 

' ' 'I- z)--^ 2Fi{-~ip,- + ip,--,l- z) 



' V{?> + ip)V{?>-ipy ' ^ ''2 "-'2'""' 2' 

Notice that for the eigenvalues of the Laplacian on S^, i.e. p = in {n > 3), the second term 
on the right-hand side vanishes. In this case the two choices are simply related by (— l)""*"^ 
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and they are both regular for all /i. Since F{1 — — ^ 1 for coincident points, we must take 
this solution in (pO|). But when we express the discrete sum (^) as a contour integral, to 
maintain regularity of the integrand at /x = vr we need to first replace F{1 — z) by a term 
l)""*"^, and then continue the latter term to — (coshpvr)"^ 2-^1 (3 + ip,3 — ip 



Now we write the sum in (|20| ) as an integral along a contour Ci encircling the points 
p = 3i,4i, ...Ni on the imaginary p-axis, where tends to infinity. Using the analj^ic 
properties of the terms in the discrete sum extended into the complex p-plane we have for 
X > X\ 



dp gp{X)^ljp{ X') 
Ci psinhpiT a_p 



a' a! 



or ' (/i, X, X') = ^ i ^ ^(J ' (/^)- (23) 



where W^(p)'*"' (z^) is defined in the Appendix, equations (A6), (AlO) and (All), using the 
forms regular at /i = vr i.e. a{z) and I3{z). From the explicit forms given, it is clear that 
^(pT analytic in complex p-plane in the required region. To see that ( p3D is equivalent 

to the sum ( ^Of ) introduce 1 = coshpvr/ coshpvr into the integral. Then note that cothpvr has 
residue vr"^ at every integer multiple of i. Finally, use (|23|) at p = ni to rewrite (a*) 
in the form regular at /i = 0. The factor of coshpvr from cancels that in the integrand. 
Note also the minus sign that appears in (p3|) cancels that introduced by the change in sign 
of the normalisation factor Qp = —{jP' + 4)/(307r^), which is positive ii p = in, n > 2, but 
negative for real p. The cancellation of these signs ensures that the Lorentzian correlator 
has the correct positivity properties. 

We now distort the contour for the p integral to run along the real p axis (Figure 2). At 
large imaginary p the integrand decays exponentially and the contribution vanishes in the 
limit of large N. However as we deform the contour towards the real axis we encounter two 
poles in the cothpvr factor, the latter ai p = i becoming a double pole due to the simple zero 
of a_p. For the p = 2i pole, we note that it follows directly from the the normalisation factor 
Qp that Vr(2i)'' ~ 0. Indirectly, this is a consequence of the fact that spin-2 perturbations do 
not have a monopole or dipole component. At p = i we have a double pole. However, the 
bound state wavefunction is just proportional to b{X) and the metric tensor perturbation 
tij = b~^{X)tij is therefore independent of X. The latter coordinate continues to conformal 
time in the open universe, and it follows that the metric perturbation is time-independent 
and will not contribute to the Sachs- Wolfe formula (0). However to understand this mode 
more deeply, recall that for = — 1 a degeneracy appears between p^ = —1 tensor- type 



perturbations and p^ = —4 vector-type perturbations |Tl[]. To be more precise, the trace- 



lesss transverse tensors qf- '^^^'' may be constructed from the vector harmonics y^^^*)*''™' by- 
symmetrised covariant differentiation. One therefore has (fi^ij?' = — 1) = '^li\j){pi = ^4). 
This means that this discrete tensor mode is not invariant under (vector) gauge transfor- 
mations. It may be generated by a purely spatial gauge transformation without disturbing 
the value of the scalar field [O]. We may therefore use the remaining gauge freedom in the 
decomposition (|^) to set W]^^^-' = 0. We conclude that up to a term involving a pure gauge 
mode, we can deform the contour Ci into the contour C shown in Figure 2. Since the inte- 
grand involves a factor (psinhpvr)^^ which has a double pole at p = 0, we leave the contour 
avoiding the origin on a small semicircle in the upper half p-plane. We shall see that for the 
Coleman-De Luccia and Hawking- Turok cases the complete integrand is actually regular at 
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FIG. 2. Contour for the Correlator. 



p = 0, but for perfect de Sitter space the double pole survives. The contribution to the 
Green function from the small semicircle acts to regulate the integral /q°° dp/p^ coming from 
the real axis. Thus in our treatment, even in perfect de Sitter space the Green function 
is finite, in contradiction to the conclusion reached in treatments based on mode-by-mode 
matching. 

IV. TWO-POINT TENSOR CORRELATOR IN AN OPEN UNIVERSE 

The analytic continuation into the open universe is given by setting Vt = —ix and a = it 
(see and letting a(t) b{a) = —ia{ia). Here VL is the polar angle on the three-sphere. 
For our correlator, without loss of generality we may take one of the two points, say Vt' to 
be at the north pole of the three-sphere. Then /i = f2, and /i continues to —ix- We then 
obtain the correlator in the open universe where one point has been chosen as the origin of 
the radial coordinate x- 

The background line element of the Lorentzian region is 

ds^ = -dt^ + a^{t) (dx^ + sinh^ ^^^^2) • (24) 

The conformal coordinate X continues to conformal time r as follows 

["■"^ da m , . 

Jit b{a) 2 ^ ^ 

where the conformal time r is defined via 

/ r^n da /■* dt' \ , , 

^^rl^))- '''' 

We now wish to make the substitutions fi = —ix, where x is the comoving separation on H^, 
in the open universe, and X = —i^ — r. The first continuation may be done immediately. We 
use the explicit formula for the bitensor regular at fi = n, given in the Appendix, equations 
(A6), (AlO) and (All) to write the in following p- integral for the Euclidean Green function: 

G^'^'{f^, X, X') = K — ^ L{X)g_,{X') + ^g,{X)g,{X')\ W%,^,M (27) 
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where we have used the formula (jl^) to re-express ipp terms of the Jost functions 
gp{X). The obstacle to setting X = — r — y is that the integrand of (^) contains a 
term gp{X)gp{X') ~ jf simply make the substitution X = — z| — r this would 

produce a term going as e^'^. But the bitensor defined in (AlO) and (All) involves terms 
which behave as e+*'^'^^*^\ and the two factors of e^'^ would lead to a meaningless divergent 
integral. To circumvent the problem, we use the following identity. For X — X' > 0, we 
have 

dp9p{X)MXX^..F^p^ = Q (28) 



c p a- 



p 



where F{p) are the j9-dependent coefficients occurring in the final (Lorentzian) form for the 
bitensor given in (A12). This identity follows from the analyticity properties of the integrand 
explained above, and the fact that, despite first appearances, the formulae ^412 are actually 
analytic ai p = i. We now insert 1 = sinhpvr/ sinhpvr under the integral, to show that 
the integral (p8|) with a factor e^'^ / sinhpvr inserted equals that with a factor e~^'" / sinhpvr 
inserted. The resulting identity allows us to replace the dangerous terms in the bitensor 

We now perform the X continuation. The analytic continuation of the Euclidean mode 
functions is given by 

g^p{X) ^ e^gi^ir) (29) 

where the Lorentzian Jost function g^ij) is the solution to the Lorentzian perturbation 
equation Kgp^r) = p^gp^r) obeying (7^ (r) e"*^"^ as r — > — oo. Equation (p9| ) follows 
by matching at large X. We finally obtain the Lorentzian tensor Feynman (time-ordered) 
correlator, for r' — r, 

Gl-'^'ix, r, r') = . I {e-'^-giMp{r') + ^^9^ (r) g^ {r'^j W%,ix) (30) 

where the Lorentzian bitensor H^^]f/j/(x) of relevance in the hyperbolic universe is defined 
in the Appendix, equation (A12). The factor a^ja^p is simply a phase, since for real p the 
Euclidean wavefunction is real so a* = a_p. 

Now we would like to represent the result (^) as an integral over real < p < oo. 
The term cothpvr in the integrand seems to produce a double pole at p = 0. However, 
for either the Coleman-de Luccia or Hawking- Turok instantons, the refiection term in (|30D 
turns out to precisely cancel the first term as p — » 0. This cancellation seems to have first 



been discovered in refs. [g. The reason for the cancellation is that for any potential 
except a perfectly refiectionless one, at very low momenta (i.e. very long wavelengths) the 
wavefunction is completely refiected. This means that in the small p limit both ap/a_p and 
Cp/dp tend to minus one [^. By analyticity, we expect them to vanish as p^, which makes 
the integrand of (^) analytic as p ^ 0. It is however clear from the form of the potentials 
(Figure 1) that the Coleman-De Luccia instantons are much closer to the perfect non- 
refiecting solution. Therefore we may expect the regime Cp/dp —>■ —1 to set in at much lower 
p than in the Hawking- Turok case. This will lead to a larger contribution to the large angle 
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microwave anisotropies. As mentioned above, a virtue of our treatment seems to be that 
even the de Sitter result is finite. 

In the cases of interest therefore there is no singularity at p = 0, and we may take the 
contour to run along the real p-axis. Using the symmetry p —p, The right hand side of 
(Wi) becomes 



^ r -W^-]r'kx)fcothp7r[^,^(r)^^,(r')+^?",(r)^^^(r') 

Z J — oo P \ 



+ 



1 



sinhpvr 



Cl—p dp 



■ (31) 



For real p, g^p^r)^ is the complex conjugate of gpi^r) and a_p of ap. So the second term is 
imaginary but the first and third terms are real. In fact it is straightforward to see that if 



we apply the Lorentzian version of the perturbation operator K to (|3T|) with an appropriate 
heaviside function of r — r', the imaginary term will produce the Wronskian of g^p{T)^ and 
5fp(r), which is proportional to p, times 6{t — r'). Then the integral over p produces a 
spatial delta function. From this one sees that our Feynman correlator obeys the correct 
second order partial differential equation, with a delta function source. The delta function 
goes from being real in the Euclidean region to imaginary in the Lorentzian region because 
of the factor ^fj in (pT]). 

For cosmological applications, we are usually interested in the expectation of some quan- 
tity squared, like the microwave background multipole moments. For this purpose, all that 
matters is the symmetrised correlator {{tij{x),ti>j>{x')}) which is just the real part of the 
Feynman correlator. It also represents the 'classical' piece, which in the situations of interest, 
where occupation numbers of modes are large, is much larger than the quantum piece. 

For the tensor correlator we also need to restore the factor ^"^(t) to tij. It is convenient 
to define the eigenmodes $p(t) = g^{T) / a{T). The symmetrised correlator is then given by 

({t.,(x),t.,,(x')}) = 2n^j^ ^ (cothp.$^(r)$-_,(r') + ^^^^^^) W%,{x). (32) 



where W^^^J,,{x) is defined in the Appendix, equations (A3) and (A12). 



In this integral the bitensor W^^^^ (fi) equals the sum of the rank-two tensor eigenmodes 
with eigenvalue Ap = — (p^ + 3) of the Laplacian on H^. At large p, its coefficient functions 
wj^'' (see Appendix A) behave like p sin pfi. Hence the above integral converges at large p 
for both timelike and spacelike separations. Equation (|32|) is our final result for the tensor 
spectrum from singular instantons. As in the scalar calculation , and as mentioned above, 
for Coleman-De Luccia instantons the phase ap/a_p gets replaced by Cp/dp, which is the 
refiection amplitude for waves incident from X = +oo in the Euclidean region. 

Before moving on to the observational consequences of (^) we would like to make one 
more technical comment. We mentioned already that a degeneracy appears between p^ = 
tensor modes and p^ = —4 scalar perturbations. These discrete modes were initially inter- 
preted as bubble wall fluctuations p3| , p^ . However, in our approach they do not contribute 
in the scalar calculation (for / > 2) because the corresponding spherical harmonics are sin- 
gular and overcomplete on the Euclidean three-sphere. More recently the wall fluctuations 
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were argued to have re-appeared as a long-wavelength continuum contribution on top of the 



usual continuous spectrum of even parity gravitational wave modes |T^. In this way, the 



bubble wall fluctuations were found to regularize the tensor spectrum, thought to be infrared 



divergent in pure de Sitter space ||10[. Our result for the correlator for a Coleman-De Luccia 
model is indeed infrared finite and the cancellation caused by total reflection of low momen- 
tum modes allowed us to represent the result as an integral starting at p = 0. However we 
do not agree that the presence of the bubble was needed to regularize the spectrum. In our 
method, even in perfect de Sitter space we obtain a finite result, because the contribution 
of the small semicircle on the contour C shown in Figure 2 regularises the final answer. So 
in our approach the tensor spectrum in perfect de Sitter space appears to be infrared finite, 
contrary to the findings of earlier works. This issue is academic for present purposes, but 
clearly deserves further study. 



V. IMPLICATIONS FOR THE CMB-ANISOTROPY 



Gravitational waves provide an extra source of time- dependence in the background in 
which the cosmic microwave background photons propagate. The contribution of gravita- 
tional waves to the CMB anisotropy is given by the integral in the Sachs- Wolfe formula 



1 25 



ST, 



sw , 



T 



1 

2, JTe 



X=TO-T 



(33) 



where tq and Xg are respectively the observing and last scattering time for the photons and x 
is the comoving radial coordinate. The anisotropy is characterised by the two-point angular 
correlation function C(7), where 7 is the angle between two points on the celestial sphere. 
It is customary to expand the correlation function in terms of Legendre Polynomials as 



C(7) 



^(0)^(7) 



°° 2/ -I- 1 

1=2 



(34) 



where in standard notation Ci = {\aim\'^). Hence, inserting the Sachs-Wolfe integral into 
(P^) and substituting (R^ for the two-point fluctuation correlator yields 



1 fTo rro Q Q 



(35) 



In order to obtain Ci we write the bitensor back in terms of its defining tensor eigenmodes 
on H^. Since q^^"''"^ = 0, only the even parity modes contribute to the CMB-anisotropy. 
The normalised eigenfunctions q^^^'^ix-, ^) 0) can be written as Qyy{x)Yim{d, 



where 26 



p2(j)2 _|_ 



(sinh x) 



1-2 



d 



i+i 



sinhxc^X, 



(cospx) 



and 



Ni{p) 



{l-l)l{l + l){l + 2) 

vrn;=2(j'+p') 



1/2 



(36) 



(37) 
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Hence we obtain for the power spectrum of multipole moments 



Jo 2p JTe 



TO /-TO 

dr 



dr' (cothpTT [<i>J(r)<i>^p(r') 



1 



sinh pn 



(38) 



The contribution to the multipole moments due to the second, reflection term falls ex- 
ponentially with increasing wavenumber. However in contrast with the scalar fluctuations 
the long-wavelength tensor perturbations do give a substantial contribution to the CMB 
anisotropics. Hence the dependence of the tensor spectrum on the boundary conditions for 
the perturbations defined by the instanton background - Dirichlet for Hawking- Turok, free 
boundary conditions for Coleman-De Luccia, may provide a way to observationally distin- 
guish different versions of open inflation. From the discussion above, we expect a larger 
contribution at low p for regular instantons. We shall perform the numerical computation 
of the needed reflection coefficients in future work [0 . 

In addition, for a complete calculation of the Ci one must evolve the Lorentzian mode 
functions $p(r) forward from the beginning r = — oo of inflation inside the open universe 
up to the present time r = tq. In the inflationary phase of the open universe the mode 
functions closely follow perfect de Sitter evolution in which they tend to a constant after 
the physical wavelength has been stretched outside the Hubble radius. The amplitude and 
phase of this constant defines initial conditions for the radiation and matter dominated 
eras in which the modes of interest re-enter the Hubble radius. The radiation and matter 
evolution is straightforward to study numerically, and from this one can compute the Sachs- 



Wolfe integral (BSf) and the the multipole moments Ci. 



VI. CONCLUSION 

We have computed the spectrum of tensor perturbations predicted in open inflation, 
according to Euclidean no boundary initial conditions. The Euclidean path integral unam- 
biguously specifies the tensor correlators with no additional assumptions. We feel that the 
present work places earlier results on a substantially firmer footing. Our final result for 
the correlator, (P^, and the cosmic microwave multipole moments ( P^ ) is given in terms of 
scattering amplitudes in the Euclidean region and mode functions in the Lorentzian region. 
Both are straightforward to compute numerically, and we shall do so in future work 0. 
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APPENDIX A: MAXIMALLY SYMMETRIC BITENSORS 



A maximally symmetric bitensor T is one for which a*T = for any isometry a of the 
maximally symmetric manifold. Any maximally symmetric bitensor may be expanded in 
terms of a complete set of 'fundamental' maximally symmetric bitensors with the correct 
index symmetries. For instance 



Tiji'f = ti{fi)gijgi>y + t2{fi) 



9ii'9jji + 9ji'9 



+ ti{^)ninjn^,nj, 



(Al) 



where the coefficient functions tj{fj,) depend only on the distance fi{fl, Q') along the shortest 
geodesic from Q to Q'. n^,{Q,fl') and n^{Q,Q') are unit tangent vectors to the geodesies 
joining Q and Q' and gijr{Q,Q') is the parallel propagator along the geodesic; V^gj is the 
vector at Q' obtained by parallel transport of along the geodesic from Q to Q' ||18|| . 
The bitensor 



(A2) 



■Plm 



appearing in our Green function (|TB|) has some additional properties arising from its con- 
struction in terms of the transverse and traceless tensor harmonics q. 



(v) 

The tracelessness of Wi^.u, allows one to eliminate two of the coefficient functions in 



(p)Vlm 



(or 



It may then be written as 



w 



ip) 



9i 



3ni 



9i- 



nigji'rij, + njgii,nj, + nig^^n^, + n^gi^iTii, + iniUjUimj, 
9ii'9jj' + 9ji'9ij' - '^n^g^,-,n - - 2n^,g^-n-, + Qn^n-n^,n-, 

ip) 



(A3) 



The requirement that the bitensor be transverse ^^Wlj/,j, = and the eigenvalue condition 



(A3 — Ap)W^(-p^-'* = impose additional constraints on the remaining coefficient functions 
w'f\jj). To solve these constraint equations it is convenient to introduce the new variables 




f a(/i) 



(p^+9)sin/x d/i 



where /i is the geodesic distance on S^. In terms of a new argument z 
transversality and eigenvalue conditions imply for a{z) 



z(l 



. dPaiz) 



da{z) 
dz 



(/ + 9)a{z) 



(A4) 

cos^(/i/2) the 
(A5) 



and then for the coefficient functions 
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wi = Qp [2{\y - Q)z{z - 1) - 2] a{z) + f + Q)z{z - \){z - l)\p{z) 

W2 = Qp 2(1 - z) [(Apr2 - Q)z + 3] a{z) - | [(A^r^ + 6)^(2 - 1)(^ - §)] f3{z) (A6) 
_ «;3 = [-2(Apr2 - Q)z{z - 1) + 3] a(^) - | [(A^r^ + 6)z(z - - 1)] (3{z) 

with Ap = (p^ + 3) on 5'^ and Qp a normalisation constant. 

To fix the normahsation constant Qp we contract the indices in the coincident hmit 
z ^ 1. This yields [|19 



(P) ij 



(fi, fi) = ^ gjf ^'"(n)g(P)^'™ '^{ny = 30Qpa(l). (A7) 



vim 



By integrating over the three-sphere and using the normalisation condition (0) on the tensor 
harmonics one obtains Qp = —3^3^^- 

Notice that ( [A5| ) is precisely the hypergeometric differential equation, which has a pair 
of independent solutions a{z) = 2-^1 (3 + ip,3 — ip, 7/2, z) and 

a(l — z) = 2-^1(3 + ip,3 — ip, 7/2,1 — z). The former of these solutions is singular at 2; = 1, 
i.e. for coincident points on the three-sphere, and the latter is singular for opposite points. 
The solution for f3{z) follows from (|A4|) and is given by 



P{z)= 2F,{A-tp,A + ip,9/2,z). (A8) 
The hypergeometric functions are related by the transformation formula (eq. [15.3.6] in P7| ) 

T{c)T{c-a-b) 



2Fi{a,b, c,z) 
r{a + h- 



T{c-a)T{c-h) 



Fi{a, b, a + b — c,l 



_, ^(1 - z) 2Fi{c - a,c - b,c - a - b,l - z). 



(A9) 



Only for the eigenvalues of the Laplacian on S^, i.e. p = in {n > ?>), the term on the 
second line vanishes for 2-Fi(3 + ip,?> — ip, 7/2, z). In this case the functions are related by 
(— l)"""*"^ and they are both regular for any angle on the three-sphere. But since F{l — z) —>■ 1 
for coincident points, we must choose a{l — z) in the bitensor appearing in the Eucliden 
Green function (PP]). This choice also follows from matching the delta function in the Green 
equation itself. In fact, the hypergeometric series terminates for these parameter values and 
the hypergeometric functions reduce to Gegenbauer's Polynomials C^_3(l — 2z). 

We conclude that the above properties required of the bitensor completely determine its 
form. Notice that in terms of the label p we have obtained a 'unified' functional description 
of the bitensor Wj^^ on and although its explicit form is very different in both cases. 
In fact it is precisely this which allowed us in Section IV to analytically continue the angular 
part of the Green function from the Euclidean region into the open universe. 

To perform the continuation we note that the Euclidean geodesic separation /i continues 
to —ix where x is the comoving geodesic separation on H^. We apply the relation ( |A9| ) in 
an intermediate step of the calculation, the continuation of the bitensor into the complex 
p-plane. In this step the functions a{z),(3{z) rather than a(l — z) and (3{1 — z) enter. The 
hypergeometric functions on are defined by analytic continuation (eq. 15.3.7 in [^) and 
may be expressed in terms of associated Legendre functions as 
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aiz) 



f (-sinhx)~^/^P 
f(-sinhx)"^/^P 



-5/2 , 
-1/2+tpl 
-7/2 , 
-1/2+ipl 



cosh x), 
cosh x). 



(AlO) 



Using the relation — cosh(x) = cosh(x — zvr), the Legendre functions may be expressed as 
P:i}2+ip{-coshx) = 

Pllj^+ipi-coshx) = 



3^^(1 + p2)-i(4 +p2)-i [3 cothxcoshp(7r + ix) 

isinhp(ix+TT) ( ^r, _2^n i „^+u2,A , , ^2^„„„„„^,2, 



2p 



(2 - p^){l + coth^ x) + (4 + p2)cosech''x 



'1 



p2)-i(4+p2)-i(9 + p2)-i^ 



TT sinh X 

coshp(7r + ix){p^ — 11 — IScosech^x) 

_gismhp(!X±l) hi _ p2) pQ^l^S ^ + (p2 ^ 3) ^^^j^ ^ cOSech\ 



(All) 



The factors e^^'^ in these expressions combine with similar factors from the continuation of 
the conformal spatial coordinate X to produce our final result (P^). The coefficient functions 



of the bitensor W^M.i (x) in our final result (p2[) for the tensor correlator are 



W2 



cosech^X 
47r2(p2+i) 



smpx( 
P 



47r2(p2+i) 



(3 + {p^ + 4) sinh^ x ~ p^(p^ + 1) sinh^ x) 
+ cos px (3/2 + (p^ + 1) sinh^ x) sinh 2x) 

(3-12 cosh X - 3p2 ( 1 - 2 cosh x) sinh^ x 

12 + 3 coshx 



_l_p2^p2 _|_ sinh^ -1^) — cospx( 

+2(p^ — 2) sinh^ x ~ 2(p^ + 1) cosh x sinh^ x) sinh x 



(A12) 



W3 



cosech^ X 
47r2(p2 4.i) 



(3 — 3p^ sinh^ x + P^{p'^ + 1) sinh^ x) 
— cospx(~3/2 + (p^ + 1) sinh^ x) sinh 2x) 

As mentioned before, for x ^ these functions converge and they exponentially decay at 
large geodesic distances. We also mention that in this form one should take the normalisation 
factor Qp to be positive, as explained in the text. 

Finally, let us mention that the scalar Green function ||^ may also be described in terms 
of hypergeometric functions. In terms of the variable z, the equation for the angular part 
Cp{ii) of the scalar Euclidean Green function (eq.(35) in P) reads 



zil-z) 



+ 



3z 



"^^'^'^ - ip' + l)Cp{z) 



dz 



(A13) 



If we express the Green function as an infinite sum (eq. (38) in 
regular at /i = and /i = tt is 



), the appropriate solution 



Cpiz) 



QpF{l + ip, 1 - ip, 3/2, 1 



Qp sinh pn 
psin n 



(AM) 



As for the tensor correlator, the normalisation constant Qp is determined by the normal- 
isation of the scalar harmonics on S*^. However, because of the extra factor (A3 + 3) in 
the scalar Green equation (eq.(35) in [§), we must also divide by 4 + in this case. This 
reproduces precisely the angular part of the scalar Green function (eq.(38) in |^). 

When expressing the Euclidean Green function as an integral, we continue Cp{z) into 
the complex p-plane, and again need to express it in terms of the hypergeometric function 
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regular at 2; = 0. We re-express F(l + ip,l — ip, 3/2, 1 — z) using the relation ( |A9|) and 
obtain 



sinhpu sinhpfvr — u) coslipu / , , 

cothpvr — = — + —. (A15) 

p sin /i p sinli pn sin /i p sin /i 

The factor cothpvr is needed in converting the sum into a contour integral. The first term is 
regular for opposite points and leads exactly to the angular part of the Lorentzian correlator 
(eq.(46) in 0) in the same way as described above for tensor fluctuations. The second term 
is a bit more subtle. Its analogue in the tensor correlator did not contribute to the contour 
integral because it had no poles within the contour. However, in the scalar we need to take 
into account the extra normalisation factor -^^^ which has a pole at p = 2i. This is the 
underlying reason for the presence of the extra term in the integral representation of the 
scalar Euclidean Green function (2nd term in eq.(37) in |^). As explained in |^), the {n — fl) 
factor in front of arises from matching the delta function in the Green equation (eq.(35) in 
0], which unlike the tensor Green equation is fourth order in derivatives. This is also the 
we had to include the extra factor Nevertheless it is clear that the scalar and tensor 

cases are very closely parallel. 
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